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Abstract 

Motivated by the recent measurement on B{Bs — )• <j)fi~^fi~) by CDF collaboration, we 
study the effects of a family non-universal Z' boson on rare semileptonic Bg — t- 
decay. In our evaluations, we analyze the dependences of the dimuon invariant mass 
spectrum and normalized forward-backward asymmetry on Z' couplings and show that 
these observables are highly sensitive to new Z' contributions. Three limiting scenarios are 
presented in the detailed analyses. Numerically, within the allowed ranges of Z' couplings 
under the constraints from Bg — Bg mixing, B — )• ttK, B^ — )• (Xg, K, K*)i^l^ fj,^ decays 
and so on, B{Bg — t- (pji^ and A^p^{Bg — > (pji^ could be enhanced by about 96% 
and 17% (133%) respectively at most by Z' contributions. However, B{Bg — )■ (p^'^ is 
hardly to be reduced. Furthermore, the zero crossing in Aps^Bg — )■ (f)fj,~^fi^) spectrum at 
low dimuon mass always exists. 
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1 Introduction 



Rare B decays induced by the flavor-changing neutral current (FCNC) occur at loop level in the 
Standard Model (SM) and thus proceed at a low rate. They can provide useful information on 
the parameters of the SM and test its predictions. Meanwhile, they offer a valuable possibility 
of an indirect search of new physics (NP) for their sensitivity to the gauge structure and new 
contributions. Experimentally, the fruitful running of BABAR, Belle and Tevatron in the past 
decade provides a very fertile ground for testing SM and probing possible NP effects. As particle 
physics is entering the era of LHC, Bg physics has attracted much more attention. 

Recently, CDF collaboration has reported the first observation of the rare semileptonic 
Bg — > (pii^ iJi" decay and measured its branching fraction to be [1] 

B{Bs 0/i V") = [1-44 ± 0.33(stat.) ± 0.46(syst.)] x 10^^ CDF . (1) 

Theoretically, many evaluations for Bg — > (j)fi'^fi~ decay have been done within both SM and 
various NP scenarios (for example, Refs. [2, 3]). The SM prediction for B^^{Bs — )■ (pfi'^fi") (~ 
1.65xlO~^(QCDSR) [2] , for example) agrees well with CDF measurement (1.44±0.57) x 10^^ for 
large experimental error. If more exact measurement on Bg — )■ (f)fi~^fi~ is gotten by the running 
LHC-b and future super-B, the possible NP space will be strongly constrained or excluded. So, 
it is worth evaluating the effects of the possible NP, such as a family non-universal Z' boson, 
on Bg — (f)fi~^fi^ decay. 

A new family non-universal Z' boson could be naturally derived in certain string construc- 
tions [4] , Eq models [5] and so on. Searching for such an extra Z' boson is an important mission 
in the experimental programs of Tevatron [6] and LHC [7]. The general framework for non- 
universal Z' model has been developed in Ref. [8]. Within such model, FCNC in 6 — s and 
d transitions could be induced by family non-universal U{1)' gauge symmetries at tree level. 
Its effects on 6 — s transition have attracted much more attention and been widely studied. 
Interestingly, the behavior of a family non-universal Z' boson is helpful to resolve many puzzles 
in B(^u,d,s) decays, such as "vr/r puzzle" [9, 10], anomalous Bg — Bg mixing phase [11, 12] and 
mismatch in ApsiB — > K*fi'^fi~) spectrum at low region [13, 14]. 

Within a family non-universal Z' model, Bg — )■ (f)fi^fi^ decay involves b — s — Z' and 
H — fi — Z' couplings, which have been strictly bounded by the constraints from Bg — Bg mixing. 
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B — )■ 7rK^*\ pK, Bd Xs/i/i, K^*^ nn decays and so on [10, 12, 13]. So, it is worth evaluating 
the effects of a non-universal Z' boson on Bg — > decay and checking whether such settled 

values of Z' couplings are permitted by CDF measurement on B{Bs ^/i^/i"). 

Our paper is organized as follows. In Section 2, we briefly review the theoretical framework 
for h — 7- sl^l~ decay within both SM and a family non-universal Z' model. In Section 3, the 
effects of a non-universal Z' boson on Bg — )■ (pfi^fi" decay are investigated in detail. Our 
conclusions are summarized in Section 4. Appendix A and B include all of the theoretical 
input parameters. 



2 The theoretical framework for b — > 5/+/ decays 

In the SM, neglecting the doubly Cabibbo-suppressed contributions, the effective Hamiltonian 
governing semileptonic b — )■ s£^£~ transition is given by [15, 16] 

^cff = --/(VtbV* J2 . (2) 

^ 1=1 



Here we choose the operator basis given by Ref. [15], in which 
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0, = ^{dj,PLb){h''l), 0,o = -^{di,PLh){lYld). (3) 

Wilson coefficients Ci can be calculated perturbatively [17, 18, 19, 20], with the numerical 
results listed in Table 1. The effective coefficients C^^J^ , which are particular combinations of 
6*7,9 with the other Q, are defined as [15] 

^eff_47r 1 4 20 80 

Cf = -C, + Yiq^), C,t = -C,o, (4) 

as as 

in which Y[q^) denotes the matrix element of four-quark operators and given by 

4 1 4 64 

= hiq^m,){-C^ + C2 + 6C■s + 60C,)--h{q^m,){7C■s + -C^ + 76C5 + -Ce) 

-h{q\ 0) {C, + ^C, + I6C5 + y Ce) + ^C, + ^C, + . (5) 

We have neglected the long-distance contribution mainly due to J/\E' and in the decay 
chain Bg — )■ (p"^^'^ (f)l~^l~ , which could be vetoed experimentally [1]. For the recent detailed 
discussion of such resonance effects, we refer to Ref. [21]. 



Table 1: The SM Wilson coefficients at the scale yU = rrih. 



Climb) C2{mb) Csirrib) Ci{mb) 




C6(mb) 


CfK) CfK)-F(g2) 




-0.284 1.007 -0.004 -0.078 


0.000 


0.001 


-0.303 4.095 


-4.153 



Although there are quite a lot of interesting observables in semileptonic h — )■ si'^C.~ decay, we 
shall focus only on the dilepton invariant mass spectrum and the forward-backward asymmetry 
in this paper. Adopting the same convention and notation as [22], the dilepton invariant mass 
spectrum and forward-backward asymmetry for Bs — )■ (f)i~^£~ decay is given as 



dr<^ 



1 



3 s 



uis] 



Ami 
A 



4m2 
1 



\B\\X 



\C\\\ 



uis] 



uis 



8m;(s + 2m^")) + |F|"(A 



+ 8m^(s - Ami)) 



\gNx 



uis] 



+ Ami{2 + 2ml- s) 



2mi 



Re{BC*){X 



uis] 



-){l-ml-s) 



+ Re(FG'*)((A 



uis] 



g -){l-m'^-s)+AmiX) 

TT) Til 

-2-4X [Re{FH*) - Re{GH*){l - mj)] + ^sX\H\ 



(6) 



ds 



Wts^tbl 



2%5 



s u{s)^ 



X 



ReiC'^'^CtrWAi + I^Re(C7'^fl^Cit*) (VT2{1 - m^) + AiTi(l + m^) 



(7) 



with s = and 5 = 8/ m^_^ . Here the auxiliary functions A , B ,C , E , F and G, with the 
explicit expressions given in Ref. [22], are combinations of the effective Wilson coefficients in 
Eq. (4) and the Bs ^ (p transition form factors, which are calculated with light-cone QCD sum 
rule approach in Ref. [23] and given in Appendix B. From the experimental point of view, the 
normalized forward-backward asymmetry is more useful, which is defined as [22] 

d^FB 



ds 



d^FB ^ dr 

ds ds 



(8) 



A new family non-universal Z' boson could be naturally derived in many extension of SM. 
One of the possible way to get such non-universal Z' boson is to include an addition f/'(l) 
gauge symmetry, which has been formulated in detail by Langacker and Pliimacher [8]. Under 
the assumption that the couplings of right-handed quark flavors with Z' boson are diagonal, 
the Z' part of the effective Hamiltonian for h — )■ sl^l~ transition can be written as [11] 



1C 

V2 L VtbVi-g VtbVts J 

With the assumption that no significant RG running effect between Mz' and scales, 
Z' contributions could be treated as modification to wilson coefficients, i.e. CgiQ^Mw) = 
Cgf^i^Mw) + ACq{Mw) ■ As a result, Eq. (9) could also be reformulated as 

H^jfib str) = ^VttV* [AC^Og + AC^oOio] + h.c. , (10) 



V2 



with 



AC[,{M^) = Df,^ = {Bji - B^) . (11) 

B^fj and Bfl'^ denote the effective chiral Z' couplings to quarks and leptons, in which the 
off-diagonal element 5^ can contain a new weak phase and could be written as |-B4|e*''^^. 
To include Z' contributions, one just needs to make the replacements 

_ 4-7r 

Cit ^ Qo' = -C(o, (12) 

as 

in the formalisms relevant to B^ 



3 Numerical analyses and discussions 

With the relevant theoretical formulas collected in Section 2 and the input parameters summa- 
rized in the Appendix, we now proceed to present our numerical analyses and discussions. 

In Table 2, we present our theoretical predictions for integrated branching fraction and 
forward-backward asymmetry of Bs — > (f)fi~^fi^ decay. Within the SM, we again find our pre- 
diction B^'^^{Bs — )■ 4'fi~^^~) = 1.46 X 10~^ is perfectly consistent with CDF measurement 
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Table 2: Predictions for B{Bs 0/i+/i~)[xlO"^] and Afb{Bs (f)fi+ n')[xlO-^] within the 
SM and the non-universal Z' model. 





Exp. [24] 


SM 


SI 


S2 


Seen. I 


Seen. II 


Seen. Ill 


B 


1.44 ±0.57 


1.46 ±0.10 


2.47 ± 1.18 


1.40 ±0.27 


2.86 


1.26 


1.92 






0.34 ±0.04 


0.56 ±0.27 


2.61 ±0.19 


0.64 


0.28 


0.44 






0.29 ±0.02 


0.51 ±0.25 


1.26 ±0.08 


0.59 


0.26 


0.39 


Afb 




25.6 ± 1.2 


19.4 ± 10.9 


24 ±0.03 


29.9 


26.6 


8.9 


^FB 




5.7 ±0.6 


6.0 ± 7.4 


0.09 ± 0.02 


13.3 


6.9 


1.4 


A" 

^FB 




34.1 ±0.2 


22.5 ± 12.9 


-0.07 ±0.01 


35.0 


34.8 


13.1 




(1.44 ± 0.57) X 10~^. The forward-backward asymmetry for Bs — )■ (pfi^fi^ decay is eval- 
uated at ~ 25%, which hasn't be measured by the experiment. In addition, in Table 2, 
we also calculate their results B^'^ and Apj^ at both low {IGeV"^ < s < GGeV"^) and 
high {lAAGeV^ < s < 25GeV^) integration regions, which are sufficiently below and above 
the threshold for charmonium resonances J/ip,ip' respectively. The dimuon invariant mass dis- 
tribution and forward-backward asymmetry spectrum are shown in Fig. 1. As Fig. 1(b) shows, 
similar to the situation in B^ — K*fi'^fi~ decay, the zero crossing exist in Afb spectrum at 
So ~ 3 GeV^, whose position is well-determined and free from hadronic uncertainties at the 
leading order in [17, 22, 25]. In 5° — )■ K*fi~^fi~ decay, the Afb spectrum measured by Belle 
collaboration [26] indicates that there might be no zero crossing, which presents a challenge to 
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Table 3: The inputs parameters for the Z' couphngs [12, 13]. 





1541 (X 10-3) 




^^,^(xlO"^) 


D^^«(xl0-2) 


SI 


1.09 ±0.22 


-72 ±7 


-2.8 ±3.9 


-6.7 ±2.6 


S2 


2.20 ±0.15 


-82 ±4 


-1.2 ± 1.4 


-2.5 ±0.9 



the SM in low s region. If the future measurement on ApsiBs — ^ 0/^^/^ ) spectrum presents a 
similar result as the one in 5° — K*fi'^fi~ decay, it will be a significant NP signal. 

Within a family non-universal Z' model, the Z' contributions to Bg — )■ (f)fi~^fi~ decay involve 
four new Z' parameters l-B^I' ^lll^ ^^m^- Combining the constraints from Bs — Bg 
mixing, B — t- ttK^*^ and pK decays, \B^^\ and 0^ have been strictly constrained [10, 12]. After 
having included the constraints from B^ Xgfifi, Kfifi and K*fifi, as well as Bg fifi decays, 
we have also gotten the allowed ranges for 5*^^ and -D^^ in Ref. [13]. For convenience, we 
recollect their numerical results in Table 3, in which SI and S2 correspond to UTfit collabora- 
tion's two fitting results for Bg — Bg mixing [27]. Our following evaluations and discussions are 
based on these given ranges for Z' couplings. With the values of Z' parameters listed in Table 3 
as inputs, we present our predictions for the observables in the third and fourth columns of 
Table 2. 

As illustrated in Fig. 2, integrated branching fraction for Bg — )■ (j)fi~^fi~ is sensitive to the Z' 
contributions. Obviously, Bg — > (f)fi~^fi~ is enhanced by the Z' contributions with large negative 
5*^^, Dj^l^ and (f)^. Moreover, compared Fig. 2 (a,b) with (c,d), we find the effects of solution 
SI is more significant than the one of S2. So, for simplicity, we just pay our attention to the 
solution SI in the following. As Fig. 2 shows, the Z' contributions with a small negative weak 
phase (j)g are helpful to reduce B{Bs — t- (f)fi~^fi~). However, because the range 0f > —65° is 
excluded by the constraints from Bg~Bg mixing and B — t- ttK decays [10, 12], B{Bs — )■ (f)fi~^fi~) 
is hardly to be reduced so much by Z' contributions. 

In order to see the Z' effect on ApB^Bg — )■ (pji'^ji") explicitly, with Y{q^) being excluded, 
we can rewrite Re{CfCl^*) and Re{CfCl^*) in Eq. (7), which dominates ApsiBg (pp^fi') 
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Figure 2: The dependence of B{Bs — 0/-i^/U ) on and -D^^ within their allowed ranges in 
SI and S2 with different (p^ values. The black dashed line corresponds to the SM result. 




(a) 



Figure 3: The dependence of Afb{Bs — > 4>fJ''^fJ' ) on 5"^^^ and D^^^ at s = 1.5 GeV^ (a) and 
s = 15GeV2 (b) with = 1.09(xl0-^), (j)^ = -72° (SI) and the central values of the other 
theoretical input parameters. The blue planes correspond to SM results. 



in high and low s regions respectively, as 
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Re{CfC!^*) ^ Re(Cf ) Re(C,t*) + Im(AC^) Im(AC;o) , (13) 

Re{CfC1^*) ~ Re{Cf)Re{C1^*) . (14) 
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Combining Eq. (11) and Eq. (14), due to the tiny Z' contribution to Cy^, the only solution to 
enhance Ap^ in low s region is a larger negative -D^^f , which also can be found in Fig. 3(a). In 
high s region, as Fig. 3 (b) shows, ApB could be reduced significantly and enhanced a bit by 
Z' contributions. 

Based on the analyses above, in order to evaluate the exact strength of Z' effects, our 
following analyses can be divided into three limiting scenarios: 

Scenario I 

In order to get the maximum B[Bs (f)fi~^fi^), within the allowed ranges for Z' couplings listed 
in Table 3, we choose a set of extreme values 

= 1.31 X 10"3 , 0f = -79° , S^;^ = -6.7 X 10^^ , Z^^^^ = -9.3 x 10"^ Seen. I , (15) 

named Scenario I. With the central values of the other theoretical input parameters, we get 
B{Bs 0yu+/i-) = 2.86 X 10-^ which is 2.5a larger than CDF result (1.44 ± 0.57) x 10"^. 
Compared with the SM prediction 1.46 x 10~^, we find B{Bs (f)fi'^fi~) could be enhanced by 
about 96% at most by Z' contributions. 

This scenario is the most helpful solution to moderate the discrepancy for Afb{Bs — )■ 
K* between SM prediction and experimental data in low s region [13, 14]. As Fig. 3 (a) 
shows, we find Scenario I also provides the most helpful solution to enhance Apb{Bs — ?■ (f)fi~^fi~) 
in low s region. Compared with the SM results, we find A^p^{Bs — )■ (f)^~^fi~) could be enhanced 
by about 17%(133.3%) at most. However, in the high s region, the effect of Scenario I on 
Afb{Bs — > (pfJ''^ as Fig. 3 (b) shows, is not significant. 

In addition, due to the strong constraints on -D^^ from Bd — ?■ Xgfifi decay, the much larger 
value l-D^^I > 9.3 x 10"^ is forbidden [12], which means the sign of Re{C^^C^Q*) can hardly 
be fiipped by Z' contributions [13]. So, as Fig. 1 (b) shows, the zero crossing in Aps spectrum 
also exists and moves to sq ~ IGeV"^ point in this scenario. 

Scenario II 

From Fig. 2, one may find that B{Bs — )■ (f)fi'^fi~) can hardly be reduced by Z' contributions so 
much within the allowed Z' parameters' ranges. The most minimal value of B{Bs — )■ (f)fi~^fi~) 
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appears at 

\Bi\ = 1.31 X 10-3 , 0f = -65° , 5^^^ = -2 X lO'^ , DJ^^ = -4 x lO'^ Seen. II , (16) 

named Scenario II. In this scenario, compared with SM prediction, we find B{Bs — (f)fi~^fi~) 
could be reduced just by about 14% at most by Z' contributions. Due to the small Z' contri- 
butions, its effect on ApsiBg — )■ (pji'^ji^) is also tiny. 

Scenario III 

As Fig. 3 (b) shows, ApB^Bg — )■ would be reduced rapidly in high s region when 

is enlarged. So, we present a limiting scenario for the minimal Ap^lBg — )■ (j)fi~^fi~), 

= 1.31 X 10-^ , 0f = -65° , = 1.1 X 10^^ , = -9.3 x 10^^ Seen. Ill , (17) 

named Scenario III. Compared with SM prediction, A^pj^{Bs — )■ ^/i'^/i") is reduced by about 
62% (62%). However, as Fig. 3 (b) shows, in the low s region, ApB is just enhanced a bit. 
So, this scenario also leads to the minimal ApBiBg — t- (pfi^fi') ~ 8.9%, which is 65% smaller 
than SM prediction. While, in this scenario, our prediction B{Bs — )■ (f)fi~^fi~) = 1.92 x 10"^ also 
agrees with CDF measurement within Icr. So, although Scenario III presents a strange effects 
on ApB spectrum, it is not excluded by current measurement either. Moreover, different from 
Scenario I, zero crossing in ApB spectrum moves to positive side in this scenario. 

4 Conclusion 

In conclusion, motivated by recent measurement on B{Bs — ?■ (f)fi^fi~) by CDF Collaboration, 
after revisiting Bg — ?■ (j)fi~^fi~ decay within SM, we have investigated the effects of a family 
non-universal Z' boson with the given Z' couplings. Our conclusions can be summarized as: 

• Branching fraction and forward-backward asymmetry for Bg — )■ (f)fi'^fi~ decay are sensitive 
to Z' contributions. All of the Z' couplings listed in Table 3 survive under the constraint 
from B{Bs — ?■ 4'fi~^fi~) measured by CDF within errors. 

• We present three limiting scenarios: B{Bs — t- tt^K'^) and A^p^{Bs (pfi^fi^) could 
be enhanced by about 96% and 17% (133%) at most by Z' contributions (Scenario I); 
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However, B{Bs — t- it~K~^) is hardly to be reduced ( reduced by 14% at most in Scenario II) 
by Z' contributions; Moreover, in Scenario III, A^^j^(Bs (pn'^ji") reaches its minimal 
value, which is 65% (62%) lower than SM prediction. 

• The zero crossing in Apb{Bs (f)fi~^fi^) spectrum always exists in the three scenarios. 

The refined measurements for the Bs leptonic decay Bg — ?■ (f)fi~^fi~ in the upcoming LHC- 
b and proposed super-B will provide a powerful testing ground for the SM and possible NP 
scenarios. Our analyses of the Z' effects on the observables for Bg — )■ (j)fi'^fi~ decay are useful 
for probing or refuting the effects of a family non-universal Z' boson. 

Acknowledgments 

The work is supported by the National Science Foundation under contract Nos. 11075059, 
10735080 and 11005032. 

Appendix A: Theoretical input parameters 

For the CKM matrix elements, we adopt the UTfit collaboration's fitting results [28] 

p = 0.132 ±0.02 (0.135 ±0.04), r/ = 0.367 ± 0.013 (0.374 ± 0.026), 

A = 0.8095 ±0.0095 (0.804 ±0.01), A = 0.22545 ± 0.00065 (0.22535 ± 0.00065). (18) 

As for the quark masses, we take [29, 30] 

niu = irid = nig = 0, rric = 1.61^*^0 52 GeV, 

rrib = ^.79tom GeV, = 172.4 ± 1.22 GeV. (19) 

Appendix B: Transition form factors from light-cone QCD 
sum rule 

In order to calculate the Bg — decay amplitude, we have to evaluate the Bg (p matrix 
elements of quark bilinear currents. They can be expressed in terms of ten form factors, which 
depend on the momentum transfer between the Bg and the mesons {q = p — k) [23]: 
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Table 4: Fit parameters for 5^ — )■ </> transition form factors [23]. 





m 


ri 


m\ 


r2 








0.434 


1.484 


5.322 


-1.049 


39.52 


Eq. (22) 


^0 


0.474 


3.310 


5.28^ 


-2.835 


31.57 


Eq. (22) 




0.311 






0.308 


36.54 


Eq. (24) 




0.234 


-0.054 




0.288 


48.94 


Eq. (23) 


rpBs-^4' 


0.349 


1.303 


5.322 


-0.954 


38.28 


Eq. (22) 


rpBs^l/) 
-'2 


0.349 






0.349 


37.21 


Eq. (24) 


rf:B^<j) 


0.349 


0.027 




0.321 


45.56 


Eq. (23) 



{m\dl,{l-l^)h\B,{p)) = -ie*(mB,+m^)A^{q^) + i{2p-q),{e* -q) 



with A;{q') = ^^^A,iq') -^^^A^iq') and Ao(0) = ^3(0), 



(20) 



{m\s(J^.uq''{l + l,)b\Bs{p)) = 2e;..p.e*Vr2Ti(g2) 

+T2(g2) [e;(m|^ - mj) - (e* ■ q) {2p - q). 



-Uq^W-q) 



5m 



m 



B., 



{2p-q). 



(21) 



with Ti(0) = T2(0). is the polarization vector of the (f) meson. The physical range in s = 
extends from Smin = to Smax = (m^^ — 'm^Y . 

These transition form factors have been updated recently within the light-cone QCD sum 
rule approach [23]. For the q^ dependence of the form factors, they can be parameterized in 
terms of simple formulae with two or three parameters. The form factors V , Aq and Ti are 
parameterized by 



Fis) = 



ri 



s/m? 



+ 



r2 



R 



1 — s/m 



(22) 
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For the form factors A2 and T3, it is more appropriate to expand to the second order around 
the pole, yielding 

^(^) = 1 '\ 2 + TT^V^ ' (23) 

1 — s/m^ (1 — s/my 

where m = mgt for A2 and T3. The fit formula for Ai and T2 is 

9 9 

^^-^ — ftl r . ~. _ -1 

The form factor T3 can be obtained through the relation T^is) = K3(s) -^2(5) . All the 

relevant fitting parameters for these form factors are taken from Ref. [23] and are recollected 
in Table 4. 
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